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Darquenne, Chantal, Peter Brand, Joachim Heyder, 
and Manuel Paiva. Aerosol dispersion in human lung: 
comparison between numerical simulations and experiments 
for bolus tests. J.Appl. Physiol. 83(3): 966-974, 1997.-Bolus 
inhalations of 0.87-^m-diameter particles were administered 
to 10 healthy subjects, and data were compared with numeri- 
cal simulations based on a one-dimensional model of aerosol 
transport ancrdeposition in the human lung {J. AppL Physiol. 
77: 2889-2898. 1994). Aerosol boluses were inhaled at a 
constant flow rate into various volumetric lung depths up to 
1,500 ml. Parameters such as bolus half-width, mode shift, 
skewness. and deposition were used to characterize the bolus 
and to display convective mixing. The simulations described 
the experimental results reasonably well. The sensitivity of 
the simulations to different parameters was tested. Simu- 
lated half-width appeared to be insensitive to altered values 
of the deposition term, whereas it was greatly affected by 
modified values of the apparent diffusion in the alveolar zone 
of the lung. Finally, further simulations were compared in 
experiments with a fixed penetration volume and various flow 
rates. Comparison showed good agreement, which may be 
explained by the fact that half-width, mode shift, and skew- 
ness were little affected by the flow rate. 

convective mixing; bolus inhalations; aerosol transport simu- 
lation 



several authors (8, 10, 14, 16-18) have shown that, 
during respiration, fresh air is irreversibly transferred 
from the inspired air into the surrounding air. This 
transfer is attributed to Brownian diffusion and convec- 
tive mixing. Convective mixing refers to all the mecha- 
nisms, except Brownian diffusion, that are involved in 
the transfer. For example, convective mixing occurs as 
a result of dispersion processes, depending on such 
factors as velocity patterns, airway and alveolar geomr 
eiry, asymmetries between inspiratory and expiratory 
flows, nonhomogeneous ventilation of the lung, and 
cardiogenic mixing. However, the contribution of each 
mechanism remains uncertain. Convective mixing in 
the lung has been measured by using the aerosol bolus 
technique (2, 10), whereby particles are not distributed 
over the entire inhaled volume but are confined within 
a small volume (bolus) of the ins pi red* air. Inasmuch a s 
-1-^m-diameter particl es h ave very low intrin sic mo- 
ti ons, they act as a nondiffusing g^ s (D, and they may 
be used to trace convective ancfbulk processes. The 
bolus undergoes progressively more axial dispersion as 
it passes through the respiratory tract. This dispersion 
may be easily measured by comparing aerosol concen- 
tration curves vs. volume recorded at the mouth during 
inspiration and expiration. Furthermore,. depending on 
the bolus location within the inspiratory phase, it 



reaches different zones of the lung and, therefore, 
allows a probe of convective mixing at predetermined 
depths within the respiratory tract. 

We provide experimental data of aerosol bolus tests 
and compare the data with numerical simulations 
based on a one-dimensional model of aerosol transport 
and deposition within the human lung (6). Bolus inhala- 
tions to various penetration volumes were adminis- 
tered to 10 healthy subjects. Parameters such as bolus 
half-width, mode shift, skewness, and deposition are 
used to characterize the bolus and to display convective 
mixing. Numerical computations simulating the experi- 
ments are completed, and experimental and numerical 
data are compared, Additional experimental data ob- 
tained in 79 subjects by Brand et al. (3) are used in the 
comparison. The sensitivity of parameters of the nu- 
merical model is also tested. 

Glossary 

C Aerosol concentration 

d Airway diameter 

d p * Particle diameter 

D Diffusion coefficient 

D t Apparent diffusion coefficient 

D B Brownian diffusion coefficient 

DE Particle deposition (expressed in 

ex Expiratory phase 

FRC . Functional residual capacity 

g Gravitational acceleration 

H Bolus half-width 

in Inspiratory phase 

2 Airway length 

L Tbtal deposition function 

L d Deposition function due to diffusion 

Li . Deposition function due to inertial impaction 

L K Deposition function due to gravitational sedim«f n: a- 
tion 

M Bolus mode 

MS Mode shift of the bolus 

N{z) Number of airways in generation z 

iVp Number of particles 

N 9 [z ) Number of alveoli in generation z 

Q Flow rate 

RV Residual volume 

s Total airway cross section 

.s a ' Inner surface area of alveolus 

S Alveoli + airway cross section 

Sk Skewness of the bolus 

St Stokes* number (p p d\ u/18/icH 

t Time 

TLC Total Lung capacity 

u Mean axial velocity averaged over the cross section > 

u* Mean uxial velocity averaged over the cross section .< 

u„ Gravitational settling velocity 

V Volume 

r eum Cumulative volume 
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ft 
<t>d 

Hp 
a 



Penetration volume 
Axial coordinate 
Generation number 

Fraction of alveolated surface of airway - 

Deposition rate by diffusion 

Deposition rate by sedimentation 

Dynamic viscosity of air 

Center of mass of the bolus 

Particle density 

Standard deviation of the bolus 



METHODS 

Numerical simulation. The bolus experiments are simu- 
lated using a numerical model developed earlier by Dar- 
quenne and Paiva (6). In the simulations, a one-dimensional 
equation describing the aerosol transport and deposition is 
solved within a trumpet model based on the morphometric 
data of Weibel (19) and Haefeli-Bleuer and Weibel (9). This 
trumpet model is similar to model C of Darquenne and Paiva 
(6). The equation describing the aerosol transport and deposi- 
tion along the airway path is 



with 



r/C tt * 2 C 1 S{sD) SC 

— = -£> + 

in S fix 2 S i)x dx 



QaC 

S dx 



L 

S 



(2) 



where C is aerosol concentration, t is time, D is diffusion 
coefficient, s is total airway cross section, S is alveoli + airway 
cross section, x is axial coordinate, Q is flow rate, and L is a 
deposition term that, incorporates deposition due to inertia! 
impaction \L\), gravitational sedimentation (L»), and Brown- 
ian diffusion (L d ). These functions are listed in the appendix. 
D incorporates Brownian diffusion (Z?b) and convective mix- 
ing (A,) and is expressed by 



D = D# + 2.400 cnv7s 
D = D H + 0.167/z/ 



V eum > 49.3 ml 



(3) 



Four parameters are used to characterize our simulations: 
change in H, mode shift (MS), skewness (Sk), and particle 
deposition (DE). H and MS are defined as differences between 
inhalation and exhalation data, whereas Sk is calculated 
from the aerosol concentration curve during expiration. 

The change in H reflects aerosol dispersion and is ex- 
pressed as 



where H xn and i* M are H of the inspired and expired boluses, 
respectively^ The MS is defined as the difference between the 
position of the expired bolus mode (M„) and the position of 
the inspired bolus (V p ) 



MS = M ux - V„ 



(5» 



A positive value of MS indicates that the position of the mode 
of the expired bolus has shifted to a larger lung volume than 
its location in the inhaled bolus. 

Sk is computed from the expiratory curve and is expressed 

by 



Sk = 



1 /(V - £) a C(V)dV 
a * J C(V) dV 



with 



/ V C(V) dV 

£ = _ 

J C(V) dV 



[6) 



where V and C(V) are the volume and the aerosol concentra- 
tion, respectively, and with 



// (V - p) 4 C(V)dV 



(ftj 



V S C(V) dV 
Finally, the percentage of deposited particles is obtained by 



where u is mean axial velocity in the airway and I is airway 
length. The first relation is used for cumulative lung volume 
^Vum 1 <49.3 ml. i.e.. within the oral laryngeal path. The term 
2.400 cm-/s has been used previously (6) and is examined in 
the discussion. The second relation is used for V cum > 49.3 ml, 
i.e.. from the entrance of the trachea down to the alveolar 
sacs. The term 0.167u/ was proposed by Ultman (16) and 
corresponds to axial streaming for fully developed flow in a 
duct of length /. In the lung airways, the mechanisms are 
much more complex and involve a mixing process due to 
nonuniformities in the axial velocity profile in the airways. 
Also, for the particle size considered in this study, is much 
smaller than D n anywhere in the lung, except in the last 
generation, where u = 0. 

The same mesh and numerical procedure described by 
Darquenne and Paiva (6) are used in these simulations. 

Characterization of the bolus. The inhaled and exhaled 
boluses are characterized by the position of their mode (M) 
and hall-width (H) on a concentration curve vs. respired 
volume, as defined by Heyder et al. (10). M is the volume at 
which the maximum aerosol concentration occurs, and H. is 
the bolus width (in ml) between the two points of one-half the 
maximum concentration. The volumetric penetration (V p ) of 
the bolus into the lung is denned as the volume of air inhaled 
from the mode of the bolus to the end of the inspiratory phase. 



DE = 



1 - 



AT 



N ■ 



100 



(9) 



where N Ptin and represent the number of inspired and 

expired particles, respectively. The ratio N VM /N p -„, is ob- 
tained by comparing the areas of the inspired and expired 
boluses on the plot of aerosol concentration vs. volume. 

All the parameters requiring the integration of the signal 
are done between limits where C(V) si 0.15C ma ,(V), inasmuch 
as the integration is usually done with experimental tracings 
to avoid errors due to the noise of the aerosol concentration 
signal recorded during the tests. 

Breathing simulation. The first series of simulations i& 
performed with an initial lung volume [functional residua] 
capacity (FRC)] equal to 3 liters. The mouth-breathing cycle 
consists of an inspiration from FRC to 70% total lung capacity 
(TLC), then an expiration to residual volume (RV). The 
maneuver is performed at a constant flow rate of 250 ml/s. 
During the inspiratory phase, an aerosol bolus is injected at a 
preselected volume characterized by its V p from 100 to 1.500 
ml. This protocol corresponds to that performed by 10 healthy 
subjects (FRC « 3.5 ±0.7 liters) in the experimental part of 
the study. 

In the second series of simulations, V p is fixed at 600 ml for 
Q from 100 to 650 ml/s. In this scries the breathing cycle 
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consists of a Miter inspiration performed from FRC, then an 
expiration to RV. 

Experimental device. The experimental setup is similar to 
that used by Brand et aJ. (4). The relevant features consist of 
a system of pneumatic valves connected to the mouthpiece. 
The system allows selection between two inhalation channels 
(filtered air or aerosol) and an exhaust tube for expiration. By 
a computer-controlled handling of the valves, an aerosol bolus 
can be introduced at various preselected positions within a 
.clean air inspiration. The breathing flow rate is continuously 
measured with a pneumotachograph (Fleisch no. 1 tube), and 
the aerosol concentration is provided by a photometer. Fi- 
nally, the valve system, the photometer, and the pneumotacho- 
graph are heated to prevent water condensation. 

Aerosol generation. Nonhy 
tictefi with fl de 




ne ous conden _ — 

s odium chloride nuclei in a commercially a vailable - aerosol 
g enerator (Mage, Lavoro e Ambiente, Bologna, ltaryJ.""The 
aerosol is then diluted with filtered air to obtain a particle 
concentration of -20,000/cm 3 . The diameter of the particles is 
0.87 //m with a geometric standard deviation <1.15. 

RESULTS 

The subjects inhaled aerosol bolus into different lung 
depths characterized by V p from 200 to 1,500 ml. The H 
of the inhaled boluses is 20 ml. Figure 1 shows ex- 
amples of experimental and numerical bolus tracings 
at different V p . From the experimental and numerical 
tracings, we characterize the bolus by means of H, MS, 
Sk, and DE. Figure 2 displays these parameters as a 
function of V p . The experimental results obtained by 
Brand et al. (3) are also shown in Fig. 2. They result 
from bolus tests performed in 79 healthy subjects for V p 
from 50 to 800 ml. In these experiments, the breathing 
cycle was an inspiration from FRC to 70% TLC, then an 
expiration to RV. 

The experimental and numerical data showed a 
linear relationship between H and V p> reflecting an 
increase in dispersion with V p . We performed a linear 
curve fitting of our experimental data (V p = 100-1,500 
ml) and of the corresponding predictions from the 
one-dimensional model. The linear approximation was 



H = 129 cm 3 + 0.516 V D 



for the experimental data and 



H - 127 cm 3 + 0.465 V„ 



(10) 



HI) 



for the numerical predictions. 

The sensitivity of the numerical model to parameters 
such as gravity or deposition was tested. Simulations 
were performed by setting the gravitational accelera- 
tion to zero and then ignoring DE in the transport 
equation (L = 0 in Eq. 1 ). Results are shown in Fig. 2. 

Simulations are also carried out by considering sev- 
eral levels of convective mixing in the alveolar zone of 
the lung (Fig. 3). First, no convective mixing is consid- 
ered in the last four generations of the lung (D 0 = 0). 
Second, in this zone a D A equal to 10% of that used in 
the classical one-dimensional simulations is used {D a = 
0.1*0.167u/). Finally, computations are performed with 
D A = 0.5*0.167^/. These numerical data are shown in 




Volume (i) 

Fig. 1. Experimental tracings of 0.87-/zm-diaineter particles tivr 
penetration volumes (V p > = 200, 400, 800. 1.000. and 1.500 ml: flow 
rate = 250 ml/s. A: experimental data from 1 subject. B: numerical 
data obtained lrom l-dxmensional model. 



Fig. 3. The entrance of the last four generations corre- 
sponds to a V p of 544 ml. 

We checked the sensitivity of the bolus parameters r.u 
the threshold C 2: 0.15C™,. We calculated H, MS. Sk, 
and DE by using a threshold of 5, 10, and 15% of C m: „. H 
and MS were the same as those computed from cho 
volume at 0.5*C max and C max that are above the thresh- 
olds. No significant differences were found in the com- 
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Fig 2 Parameters characterizing aerosol bolus tracings as a function of V p for 0. 87 -^m -diameter particles; flow 
rale = 250 ml/6. A: half-width Of). B: mode shift (MS). C: skewness <Sk). D: deposition (DE). Symbols (means 3: SD) 
and lines refer to experimental and numerical results, respectively, a. Data from our 10 subjects; O, data from 
Brand et a). l3>; solid line, reference simulations; dotted line, simulations performed at G = 0; dashed line, 
simulations where deposition mechanisms are neglected (L « 0). 



puted DE, in contrast to Sk. The differences in Sk are 
displayed in Fig. 4, where the parameter is plotted as a 
function of V p for the three different thresholds. 

Finally, the influence of Q on aerosol dispersion is 
investigated in the second series of simulations. H t MS, 
and Sk from numerical and experimental data appear 
to be linearly related to Q. The linear regressions (y = 



rnQ + b) are displayed in. Table 1. To evaluate. Die 
weight of the slope term (/nQ) over the intercept {b), we 
computed the absolute value of the ratio mQ/6 for Q of 375 
ml/s. The computed data are displayed in Table 1.. For 
experimentally derived if, MS, and Sk, mean values are 
used in the linear regressions. The standard deviations of 
the experimental data are listed in Table 1. DE decreases 
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2GO0 




Fig. 3. Parameters characterizing aerosol bolus tracings as a function of V p for O.S7-/fm-diametcr particles: flow 
rate = 250 ml/s. Symbols l means ± SDt and lines refer to experimental and numerical results, respectively, a. Data 
from our 10 subjects; O. data from Brand et aj. (3); solid line, reference simulations; dotted line, simulations performed with 
apparent diffusion coefficient iD :t ) - 0; dashed line, simulations with D a = 0. 1«0.167iW; dot-dashed line, simulations 
with D a = 0 S*0.167«/. Entrance of last 4 generations corresponds to a penetration volume of 544 ml. 



with increasing Q. Results from numerical simulations as 
well as experimental data obtained from' 10 subjects by 
Brand et al. (3) for the same protocol are shown in Fig. 5. 



DISCUSSION 



This study concentrates on the comparison between 
numerical and experimental data of aerosol dispersion 



and deposition at various volumetric depths within the 
lung as well as the sensitivity of the simulated experi- 
ments to the different model parameters. Experiment^ 
data result from bolus inhalations administered to 10 
healthy subjects. Additional experimental data of Brand 
et al. (3) are also considered. The numerical results are 
derived from one-dimensional simulations performed 
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Fig. 4. Comparison between numerical Sk data (Sk nuin ) computed for 
different threshold of aerosol concentration (C). Dot-dashed line, C £ 
0.u5C maA : solid line, C £ O.lOC^; dotted line, C a 0.15C raft4 . 

within a trumpet model of the human lung, in which a 
one-dimensional equation describing aerosol transport 
and deposition is solved (6). As suggested previously 
(6), in our simulations we have considered D a = 2,400 
cm 2 /s to describe convective mixing in the oral- 
laryngeal path. This coefficient was chosen such that 
the numerical predictions fit the experimental data 
obtained by Anderson et al. (2) in 11 healthy subjects 
ibr V p = 100-700 ml. However, we have checked the 
validity of this coefficient by comparing numerical and 
experimental tracings of bolus from three different 
subjects inhaling at. V p = 40 ml, i.e., in the oral 
laryngeal path (Fig. 6). The comparison shows that the 
use of such a dispersion coefficient in the oral-laryngeal 
region is acceptable. 

Figure 1 illustrates the spreading of exhaled boluses 
ibr different values of V p . Numerical tracings appear to 
approximate well the experimental values. The bolus is 

Table 1. Linear regression of H, MS, and Sk on Q 



y 


v = mQ + h 






H num 


m = -7.97 x 1(J~ 3 
fc = 396.1 


7.5 x 10- :< 






m -3.29 * 10 " s 
6^402.5 


3.1 x 10- 1 


88.8 


MS num 


m« 1.67 x Hi-* 
b= -20.68 


0.30 




MS WP 


m «9.60 x JO" 2 


0.83 


41.87 


/>= -43.44 








m = -9.86 x 10' 6 


1.5X10" 2 




f>- 0.250 








m = 6.24 x Hr c 


0.16 


0.075 


6 = 0.151 







H, half-width Imb; MS, mode shift (ml); -Q, flow rate (ml/s); num. 
numerically determined: exp. experimentally determined. 



Ui 

Q 




Fig. 5. DE as a function of flow rate (Q; at \\ = 600 nil and to: 
0.87-//m-diaraeter particles. Symbols (means r SD). experimental 
data of Brand et al. (3 ); solid lines, numerical results. 

more and more dispersed as it penetrates deeper into 
the lung. During the inspiratory phase, the bulus di- 
vides into several segments that become more numer- 
ous as the bolus penetrates deeper into the lung. The 
segments recombine during expiration in such a way 
that the expired bolus is spread over a larger volume 
than the inspired bolus. This means that particles are 
transferred between the bolus and the surrounding air 
during a breath. For — 1 -/ma -diameter particles, intrin- 
sic motions are very low: the diffusion coefficient is 
0.3 X 10~ 6 cm 2 /s, and the settling velocity resulting 



0.50 



0.30 - 
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0.00 1 




200.0 



80G.0 



400.0 
Expired volume (ml) 

Fig. 6. Comparison between numerical and experimental tracings nf 
expired boluses. Solid, dotted, and dot-dashed lines, experimental 
tracings from 3 subjects; bold dashed line, numerical tracing. Partidi' 
diameter « 0.87 ;/m, Q = 250 ml/s, V p = 40 ml. 
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from gravity is —32/im/s. In a study of convective and 
diffusive gas transport in canine intrapulmonary air- 
ways, Schulz et al. (12) performed aerosol bolus inhala- 
tions with 0.86- and 2. 38-/zra -diameter particles. They 
found no systematic differences between the H of the 
smaller and the larger particles, but the settling veloc- 
ity was 7.7 times larger for 2.38- than for 0.86-/*m- 
diameter particles. The dispersion of the bolus must 
therefore be attributed to mechanisms (convective mix- 
ing) other than the intrinsic particle motions. 

The different parameters chosen to characterize the 
bolus behavior are plotted in Fig. 2 as a function of V p . 
Numerical and experimental H increase continuously 
with V p , indicating that each element of the lung 
contributes to convective mixing. Numerical values of 
H fit experimental data well. The. linear regressions 
between H and V p (Eqs. 10 and 11) showed a very 
similar intercept but a reduced slope for the numerical 
predictions compared with the experimental data. We 
also compared the linear regressions with those found 
by Heyder et al. ( 10) in a previous study performed on 
17 healthy subjects. Their protocol was such that at end 
inspiration the subjects were at 62 or 88% of their TLC. 
They found 

H = 179 cm 3 + 0.539 V p (12) 

at the higher lung volume (88% TLG) and 

H = 156 cm 3 + 0.514 V p (13) 

at the lower lung volume (62% TLC). Our tests were 
performed at 70% TLC, and the slope of the linear 
regression does not show a significant difference from 
data of Heyder et al. The intercept is, however, lower in 
our experiments than in their study As in our experi- 
ments, the data of Heyder et al. suggest a larger slope 
in the regression curve than that predicted from the 
simulations. We performed additional simulations 
where we modified the coefficient D tt in the alveolar 
zone of the lung, where its applicability is uncertain. 
Instead of using the function D a = 0.167u/, where u is 
the axial velocity averaged over the total cross section S 
of the duct, we considered the function D n = kit*, where 
is the axial velocity averaged over the cross section 
of the airways .«? (without the alveoli) and k is a 
parameter that was adjusted to obtain the best fit with 
the experimental data. The use of instead of u was 
based on the two-dimensional simulations performed 
by Darquenne and Paiva (7) showing that the flow in 
the alveolar zone of the lung is mainly confined in the 
lumen of the airways. The best fit was obtained for k = 
0.005 and led to the following regression curve 

H = 100 cm 3 + 0.542 V p (14) 

The slope of this regression curve agrees better than 
Eq. 11 with the slopes derived from the experiments 
{Eqs. 10, 22, and 13), contrary to the intercept, which is 
lower than previously reported. 

Figure 2B displays MS. For V p > 100 ml, experimen- 
tal and numerical data show negative values, meaning 
that the bolus mode is shifted to a smaller lung volume 



than its location in the inspired air. This effect becomes 
more pronounced with increasing V pj especially for tho. 
experimental results. Sk is displayed in Fig. 2C. Nu- 
merical and experimental data display positive value:?, 
indicating an extended tail toward the end of expia- 
tion. Moreover, Sk appears to be maximal at small V p 
and tends to stabilize for larger V p . Also, even if the 
curve seems to be qualitatively correct (Fig. 4), we 
probably underestimated Sk by using C ^ 0.15C max; 
inasmuch as we excluded the contribution of a long taL. 
Such a high threshold was, however, chosen to mini- 
mize the parameter's dependency on the noise level of 
the experimental signals. 

The dependence of Sk on V p might be explained by 
the effect of the convective transport vs. the diffusive 
transport of the aerosol in the respiratory tract. In the 
first generations of the lung, aerosol is mainly trans- 
ported by convection. During inspiration the velocity 
profile across the section of the airways is not uniform 
(11) and allows particles in the center of the ducts to 
reach more distal generations of the respiratory tract, 
than particles located near the walls of the airway?, 
where velocity is smaller. This causes a Sk of the bolus 
toward the distal part of the lung. During expiration, 
the velocity profile is more blunted (11), and particles in 
the center of the airways travel at a slower rate than 
during inspiration, whereas particles near the walls 
travel faster, preventing the bolus to recover its origi- 
nal shape. At larger V p , this effect is attenuated by trm 
diffusive transport, which is no longer negligible, ami 
the exhaled bolus is more symmetrical. 

Experimental and numerical deposition are corn* 
pared in Fig. 2D. Deposition increases with V p because 
of a longer residence time of the particles within tin; 
lung at deep V p . Moreover, with increasing V pj particles 
reach air spaces with decreasing dimensions, enhanc- 
ing the probability of deposition by sedimentation and 
diffusion. 

Further simulations have been performed with al- 
tered values of gravitational acceleration, deposition 
processes, or convective mixing. As shown in Fie. 2. 
simulated H and Sk appear to be insensitive t.o the 
removal of the gravitational force (g = U) and ihc 
absence of deposition (L = 0), in contrast to MS 
Interestingly, these results show that deposition, as 
simulated by the model, does not affect bolus disper- 
sion. Tbtal deposition seems therefore to be a poor 
parameter to describe the mechanisms of aerosol trans- 
port in the respiratory tract, as has been shown for u 
full inhalation of aerosols in previous one-dimensiona: 
simulations (6). On the other hand, simulated MS 
approaches zero when g .= 0 or L = 0. This suggests 
that deposition is a determining factor in the shift, or* 
the bolus mode. This is in agreement with the results o:' 
Brown et al. (5), who performed an experimental study 
on dispersion of aerosol boluses in the human lung. 
They found a. significant correlation between MS and 
DE: a mode shift toward the mouth (i.e., MS < 0) was 
associated with an increase in DE. As DE increases 
with increasing V p , MS becomes more negative for 
deeper V p . The removal of the gravitational force lends 
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to a significant decrease in DE, as expected. For 
0.87-/im-diameter particles, deposition by inertial im- 
paction is negligible, and tbe dotted line in Fig. 2D 
reflects, therefore, deposition by Browhian diffusion. 
The limitations of the one-dimensional model should, 
however, reinforce the precautions in the interpreta- 
tion of these simulations. It is indeed surprising that H 
and Sk are very insensitive to gravitational accelera- 
tion and deposition. However, further speculation on 
these comparisons should await the simulations of 
these curves with multidimensional models. 

Figure 3 illustrates that convective mixing in the 
alveolar zone of the lung largely affects all the param- 
eters except deposition. The sensitivity of altered diffu- 
sion coefficient D is marked for V p > 400 ml, i.e., for V p 
such that the bolus may enter the alveolar zone. 
Darquenne and Paiva (7) studied aerosol dispersion 
within a two-dimensional model representative of the 
alveolar zone of the human lung. They discussed the 
validity of using the classical dispersion coefficient D a 
(Eg. 3) in the one-dimensional transport equation (Eq. 
1 ) to describe convective mixing in the alveolar zone of 
the lung, whereas this coefficient is based on studies 
performed in the first generations of the bronchial tree 
(16). Their simulations suggest that this coefficient 
may probably not be directly extended to the distal 
alveolar ducts, where its use overestimates mixing. 
Their results refer to a rather small subunit of the 
acinus that is ventilated synchronously. They did not 
take into account the effect of delays or asynchrony 
induced by ventilation non uniformities between groups 
of acini or large lung regions, nor did they consider the 
chaotic mixing of flow induced by the expansion and 
contraction of alveolated ducts during breath (15). 
These effects would increase dispersion. The compari- 
son between the linear regressions of the experiments 
(Eg. 20) and the numerical predictions (Eq. 11) suggest 
that the coefficient D a we used in the one-dimensional 
model slightly underestimates convective mixing in the 
alveolar zone, inasmuch as we obtained a smaller slope 
in the linear regression of the numerical data. H is very 
sensitive, to D 61 as shown by the different simulations 
displayed in Fig. 3. 

Another interesting observation in Fig. 3 is that total 
deposition does not change significantly at any level of 
dispersion. Despite the fact that this parameter is too 
weak to allow insight into the aerosol behavior in the 
lung, Fig. 2D has shown that DE is sensitive tog. This 
suggests that experiments in hypergravity (performed 
in centrifuges) or in microgravity conditions, e.g., in 
parabolic flights or in space, can indeed shed light on 
the mechanisms of particle transport and deposition in 
the human lung. 

The influence of flow rate on aerosol dispersion is 
displayed in Table 1 and Fig. 5. H appears to be little 
influenced by the flow rate, and comparison between 
numerical and experimental H shows good agreement. 
Inasmuch as experiments are performed for constant 
inspired and expired volumes, varying the flow rate 
implies a variation of the resident time of the aerosols 
in the respiratory tract. Convective mixing, which is 
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primarily responsible for aerosol dispersion, is simu- 
lated in our numerical approach by a D a proportional to 
the mean velocity of the gas in the airway. Further- 
more, the mean displacement of particles due to diffu- 
sion may be expressed by 

Ax = \/2D n Af (75* 

where At is the mean resident time. Inasmuch as D a is 
proportional to Q and At is inversely proportional u> Q. 
Ax remains constant and aerosol dispersion (i.e. : H ) is 
not influenced by Q. 

MS and Sk are also displayed in Table 1. Except for 
the experimental MS, numerical and experimental 
data appear to be little affected by Q. If we exclude MS tf- 
at low Q «200 ml/s), the regression parameters he- 
come m = 1.2 x 10" 2 and b - -2.13, the slope of the 
regression line approaching that of MS nunv Finally, 
deposition decreases with increasing Q, as shown on 
Fig. 5. The main deposition mechanism for 0.87-j/m- 
diameter particles is gravitational sedimentation. Depo- 
sition increases therefore with increasing residence 
time, i.e., decreases with increasing Q, as the respired 
volume is kept constant. 

This is the third article in a series dealing with the 
modeling of aerosol transport and deposition in the 
human lung. The first article (6) dealt with one- 
dimensional simulations and showed that total deposi- 
tion was a poor parameter to describe the aerosol 
behavior in the respiratory tract, although the model 
satisfactorily simulated the experimental data. One ol 
the main arguments supporting this observation was 
the quasi-independence of total deposition on the level 
of convective mixing introduced in the model. The 
second article (7) concerned two- and three-dlmun- 
sional simulations of particle transport in the alvuoiar 
zone of the lung and showed that the presence of ihc 
radial alveolar septa was a major factor in the penetra- 
tion of the particles in the very periphery of the lung 
and that very large particle concentration inhomopene- 
ities are expected within any acinar duct between i h* 
lumen and the adjacent alveoli, at any moment of the 
respiratory cycle. The very different convective veloci- 
ties in the center of the duct with respect to the 
adjacent alveoli were the main reason why one- 
dimensional models assuming uniform concentration 
and velocity over the cross section of the alveolar ducts 
cannot evaluate accurately the location of particle 
deposition. Here we have shown that the one-dimen- 
sional model simulates satisfactorily aerosol disper- 
sion, which appeared to be very insensitive to gravita- 
tional acceleration and deposition. Therefore, the modei 
as developed so far seems to be suitable to com putt- 
total deposition and dispersion of the aerosol but noi. to 
locate the sites of deposition along the respiratory tract . 

In conclusion, bolus inhalations have been performed 
on 10 healthy subjects for various V p . Parameters such 
as H, MS, Sk, and DE have been used to characterize 
the bolus and to display convective mixing. Numerical 
computations based on a one-dimensional model ot 
aerosol transport and deposition in the human lung 



974 



NUMERICAL AND EXPERIMENTAL DATA OF BOLUS INHALATIONS 



have also been completed to simulate experimental 
tests. Even though a quite simplified approach has 
been used, the computations appear to describe the 
experimental results reasonably well. Numerical and 
experimental data show that irreversible processes 
occur in the bronchial tree, bringing about aerosol 
dispersion. This irreversibility of convective flow may 
be attributable to several factors. One factor may be 
differences in inspiratory and expiratory velocity gradi- 
ents: during inspiration the flow divides at each bifurca- 
tion, whereas during expiration the flow continuously 
recombines. Nonreversible secondary flows appear at 
the airway bifurcations and may increase mixing. 
Asynchrony in the expansion and contraction of "the 
acinar airways and also of larger ventilatory units may 
affect aerosol dispersion. Cardiogenic oscillations may 
also be responsible for mixing. Finally, additional nu- 
merical; data have been compared with experimental 
tests performed with 10 subjects for a fixed V p and 
various Q. Except for deposition, all the parameters 
used to describe the tests appeared to be little affected 
by Q. 



APPENDIX 

Deposition functions. The deposition term in the transport 
equation (Eq. 1 ) is 



L = L { + L, + L d 



CAi) 



where L u L if and L 4 are deposition functions due to inertial 
impaction, gravitational sedimentation, and Brownian diffu- 
sion, respectively. These functions are described by Dar> 
quenne and Paiva (6). The function describing deposition by 
inertia! impactiun is based on experimental data measured in 
a cast of the upper airways, whereas the functions of deposi- 
tion by gravitational sedimentation and Brownian diffusion 
are derived from a theoretical approach. The functions are 
expressed per unit time and unit length of the airway by 

CQ 

L, = i.3 ~ (St - 0.0001) St > 0.0001 

1 [A2) 



L. « 0 



St < 0.0001 



t'ur inertial impaction, where St is the Stokes' number 

» Pl , djju/l8/«x).by 

lN{z) v t ls 
V -ns Q 



1 — exp 



N a U) 
+ <fr. — {A3) 



for gravitational sedimentation, and by 



CQ 

L tl - — (I * a) (I - exp 



I 



36D Q LN(z) 



Q 



NJtz) 
I 



(A4) 



for diffusion, where a is fraction of alveolated surface of 
airway, N{z) is number of airways in generation 2, u % is 
gravitational settling velocity, is deposition rate by sedimen- 
tation, iV tt U) is number of alveoli in generation z t <J> d is 
deposition rate by diffusion, and s a is inner surface of alveolus. 
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